arXiv: 1502.02556v2 [hep-ph] 10 Feb 2015 


Prepared for submission to JHEP 


HIP-2014-02/TH 


On the infrared behavior of the shear spectral 
function in hot Yang-Mills theory 


Aleksi Vuorinen^ and Yan Zhu^ 

^Department of Physics and Helsinki Institute of Physics, P.O.Box 64, FI-OOOI 4 University of 
Helsinki, Finland 

Departamento de Fisica de Particulas and IGFAE, Universidade de Santiago de Compostela, 
E-15706 Santiago de Compostela, Calicia, Spain 

E-mail: aleksi.vuorinen@helsinki.fi, yan.zhu@usc.es 

Abstract: We revisit the determination of the two-loop spectral function in the shear 
channel of hot Yang-Mills theory. Correcting a technical error in an earlier computation 
and extending the result with a leading order Hard Thermal Loop resummation is seen 
to improve the infrared behavior of the quantity significantly. This makes it possible to 
straightforwardly use the result in the corresponding imaginary time correlator and the 
shear sum rule. 

Keywords: Quark Gluon Plasma, Thermal Field Theory, NLO Computations 



Contents 


1 Introduction 1 

2 Unresummed calculation 3 

3 HTL resummation 6 

3.1 Naive calculation 7 

3.2 Resummed calculation 8 

4 Results 9 

4.1 Unresummed spectral function 9 

4.2 HTL resummed spectral function 10 

4.3 Sum rule and imaginary time correlator 12 

5 Conclusions 13 

A Master integrals in the unresummed calculation 14 

A.l Definitions 14 

A.2 Evaluation of the masters 16 

A.2.1 New masters 16 

A.2.2 Discrepancy with previous results 19 

B Master integrals in the naive HTL calculation 20 

C Resummed HTL calculation 23 


1 Introduction 


The shear viscosity of the quark gluon plasma (QGP) has been identihed as one of the key 
parameters describing the medium, having particular impact on the hydrodynamic behavior 
of the matter produced in heavy ion collisions at RHIC and at the LHC (see e.g. [1-5]). 
Despite the strong experimental and phenomenological motivation, a nonperturbative first 
principles tool to predict its value is, however, still lacking, even though an extensive 
amount of work has been devoted to the topic in the weak coupling [6, 7], lattice [8-10] 
and gauge/gravity frontiers [11-13]. The main issue preventing a straightforward lattice 
determination of the parameter is its inherently Minkowskian nature: According to the 
Kubo formulae, the viscosity is available as the zero frequency limit of the corresponding 
spectral function. 


V = 


lim 

CD—)-0 


UJ 


( 1 . 1 ) 
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obtained from the imaginary part of a retarded (Minkowskian) Green’s function. 

One promising attempt to overcome the difficulties involved in the determination of 
transport coefficients is the analytic continuation of lattice results for Euclidean imaginary 
time correlators, 



do; 

—p(w)- 

TT 


cosh 


— T UJ 


sinh ^ 


0 < T < id , 


( 1 . 2 ) 


proposed and tested in [14-16]. An integral part of this program is the analytic subtraction 
of short-distance divergences from the results, i.e. obtaining as much analytic information of 
the ultraviolet (UV) structure of the spectral function as possible. This task, which should 
be carried out separately at zero and hnite temperature, is most conveniently addressed 
using the machinery of perturbation theory. 

The perturbative evaluation of thermal spectral functions becomes a complicated task 
beyond leading order, and it is only rather recently that progress in this direction has been 
achieved. For pure SU(A^) Yang-Mills theory, defined via the Euclidean action (note that 
we work in Euclidean metric throughout the paper) 


5e 
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with D = 4 — 2e, /3 = 1/T and 


(1.3) 






(1.4) 


a next-to-leading order (NLO) spectral function was first determined in the bulk channel 
in [17], building on the earlier work of [18-20] (see also [21]). Since then, the techniques 
developed in these papers have been further generalized to include the case of non-vanishing 
external three-momenta [22], motivated by applications beyond QCD, most importantly 
studies of the leptogenesis scenario. 

In the context of QGP physics, an obvious goal is to extend the bulk calculations to 
the technically significantly more tedious shear channel. This challenge has indeed been 
addressed hrst on the level of the Operator Product Expansion (OPE) in [23], and later 
by considering the full NLO spectral function in [24] (see also the related T = 0 work of 
[25, 26]). While consistent with known sum rules [27, 28] as well as the arguments of [29] 
concerning the UV behavior of various Green’s functions, the latter of these calculations 
delivered somewhat surprising results. Most importantly, it was observed that in the small- 
frequency limit, the perturbative spectral function tends to a constant. While not directly 
alarming — after all, it is known that to reach the true w —?■ 0 limit of the quantity, one 
needs to perform an elaborate resummation [7] — this implied that it was not possible to 
straightforwardly apply the result to the determination of the imaginary time correlator of 
eq. (1.2). 

The situation described above clearly calls for a more detailed study of the infrared (IR) 
behavior of the shear spectral function, utilizing a Hard Thermal Loop (HTL) resummation 
to extend the region of validity of the results of [24] to frequencies of order gT. This is not 
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only important for the sake of academic interest, but also to aid the eventual extraction of 
the shear viscosity from Euclidean lattice data and to facilitate more accurate comparisons 
with recent lattice and AdS/CFT calculations [30-36]. In addition to the HTL exercise, 
we will, however, also revisit the unresummed calculation of [24], performing a completely 
independent evaluation of the shear spectral function to test the correctness of our earlier 
results. 

As the setup of our unresummed calculation is in practice identical to that of [24], we 
refrain from presenting a lengthy introduction to the technical machinery involved in this 
part of the work. Instead, we will simply walk the reader through the necessary notations 
and dehnitions in section 2. After this, we explain the details of the HTL resummation in 
section 3, and subsequently present and analyze our results in section 4. Section 5 is finally 
devoted to drawing conclusions, while appendices A-C contain some lengthy definitions 
and technical details concerning the evaluation of the master sum-integrals encountered. 


2 Unresummed calculation 


We are interested in connected Green’s functions of specific components of the energy- 
momentum tensor of SU(Al) Yang-Mills theory 

V = , (2.1) 

denoted by 


G^u^ap{.x) = (r^i.(x) r„^(o))c, 


( 2 . 2 ) 


and in particular in the associated (momentum space) spectral functions. 




- P—>■(—i[<^-l-*0+],0) 


(2.3) 


As discussed already in [23] , the most convenient way to access the momentum space shear 
correlator Gi 2 ,i 2 (.P) U dimensional regularization proceeds by introducing the projection 
operator 




D-2 


{P^A 


+ p^pPL 


(2.4) 


where P^(P) is a usual transverse projector orthogonal to the four-vectors P and U = 
(1,0). Choosing the spatial momentum p to point in the xo-i direction, this gives 


Xi_iu,ap G^y^oi0{P) = —D{D — 2){D — 3) Gi2,i2(-P) , (2.5) 


which prompts the definition 

Grj{P) ^ G^y^apiP) ( 2 . 6 ) 

and leads to the e = 0 identity 

p^(a;) = -16/912,12 (w). (2.7) 
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Figure 1. The one- and two-loop Feynman diagrams contributing to the NLO shear spectral 
function in S\J{N) Yang-Mills theory. The curly line corresponds to the gluon field, while the grey 
blob in (v) denotes the one-loop gluon self energy. 


As explained in [24], to NLO the shear spectral function obtains contributions from 
the 1- and 2-loop graphs of fig. 1, and can thus be written in the form {d-A = — 1) 

, (2.8) 

where each term corresponds to the graph with the same index. In terms of the master 
integrals defined in appendix A, these functions read (in an arbitrary covariant gauge) 

P(i)M = - 3) _ 3)p^^,(cu) -{D- 2){D - 3)p^.(a;) ,(2.9) 

P(-)(o;) = ]^D{D - 2f{D - 3)p^o + 2D{D - 2){D - 3)p^, , (2.10) 

P{m)i^) = “ 2)(L) - 3)pj-o - “ 3)(i4 - 4)^2;! , (2-11) 

P(^^\{uj) = —D[D — 2)[D — 3)pj-o + D{D — 3)[2D — 5)pj-i — 3D [D — 3) pj-2 

IQ 1 3 

- -D{p - 2){p - 3)p^o + -^D{D - 3){2D - 25)p^, + -D{D -2){D- 3)p^o 

- ^D{D - 3){2D - 7)p^i + ^D{D - 3){2D + 21)p^2 + 2{2D^ -7D + 4)p^,, 

+ (lOD^ - 47D + 52) p 5- - 4 (D^ - 2D - 2) p g, - 2{2D^ -7D + 4)p 7- 

■^h "^h "^h 

- D {2D‘^ - IID + 15) p^i , (2.12) 

^(vi)(‘^) = “^(vii)(‘^) + D{D — 2){D — 3)pj-o + D{D — 3)pj-i + 3D{D — 3)pj2 

\ J y f b b b 

+ ~ ~ ~ ~ ~ ~ ^)Pl° 

3 25 

- - 3)Pxi - - 3)%2 - 2D{D - 3)p^ 

- [D — 3')(7D — \2^pj-y + 4D(D — 3)pj6’ + 2D(D — 3)pj-r — D[D — 3)pj-i 

■^h "^h "^h \ 

+ —D^D — 2)(D — 3)pjo + 2D[D — 3)pji + D{D — 3)pj2 

4 j j j 

+ (3D2 - 16D + 12) p^3 + I {3D^ - 16D + 12) p^4 

j 2 j 

- D{D - 6)pj- 5 - 2D{D - Q)pj-p , 


(2.13) 




P{v)i^) = ~2^i^ ~ 2)^(-D - ^)Px0 - D{D - 2){D - 3)p2;i - ~ ~ ^)Pi^ 

+ —Di^D — 2)'^[D — 3)pjo + D[D — 2)[D — S)pj-i + 3[D — 2){D — 3)pj-2 

4 “^d “^d “^d 

+ 2[D — 2)‘^(D — 3)pj-3 + —D[D — 2)[D — 3)pj-o + ~ “I" ^)Pij 

~ 2^^^ ~ “ ‘^)Pio ~ 2^^^ ~ — 2){D — 3 )pj-3 

- 2{D - 2)Vx4’ + D{D - 2)pj4 -{D- 2){3D - 3)p^^ 

“I" ~ ‘^)Px^ 4(Z) — 2)pj-q’ — 2D[D — 2)pj-e + 2[D — 2)'^pj-^ 

Z "^h "^h "^h 

— D(^D — 2)pj-7 + D[D — 2)[D — 3)pj-i — D[D — 2)[D — 3)pj-2 

- ^DiD-2)iD-3)p^3. (2.14) 

We notice that the masters I*' — cancel in the sum of all diagrams, which was to be 
expected based on their absence in the results of [24]. 

The reason we have chosen to write the spectral function in a form slightly different 
from the one used in [24] is that the above formulation allows us to separate all the IR 
sensitive masters with squared propagators to (w). This part can furthermore be written 
in an alternative form using the one-loop gluon polarization tensor 


KiQ) = 


2 eab j 


1 


{D - 2)T V + 2(Q^Q. - Q" V) 


rR^ 

D-2^ {2R- Q)^{2R-Q), 


tR\Q-RY 


2 % R^{Q-RY 

the transverse and longitudinal components of which read 


(2.15) 
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D-2 


KliQ) - ^ nsg(Q)) > ^ ngg(Q). (2.16) 
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A straightforward exercise namely shows that one can write 
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„ P(P-2)(P-3) 

P{y) (‘^J - - 1 

_ DiD-2)iD-3) ^ 

D^-1 jL|“g2(Q_p)2 

where the different integrals take the forms 

^T{Q)q‘^pl 


(T»-3)(T» + 1)5iM^!^-2(T)-2) 


Q\Q-Pr 

2 ^ ^ 1) 


QHQ-P) 


2 f ’ 


Q\Q-pr 

(2.17) 


Im 


_JQ Q^iQ .P)^ J p_).(_j[(j+i0+],O) 

, D-1 

- p^o + {D- l)p^2 - 2-—p^3 


D-l 1 

D-2^R DPP2^^l~2^^i' 


(2.18) 
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nT(Q)g^ 

Lj^ QHQ-Pr 
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, 1 D-1 

— [D — l)Px2 — - Pjo — 3' 


D-2 


[Tq Q‘^{Q-Pf 


D-2 


P—>■(—i[oj+i0+],0) 


PjO + 


3L>2 _ i2,d + 10 
2{D - 2) 


‘^Pj^ 1 

(2.19) 

Pxj ’ 

(2.20) 

Pxj ■ 

(2.21) 


The virtue of this formulation is the absence of masters with squared propagators, which 
altogether avoids the need to introduce an auxiliary mass parameter to be differentiated 
upon, as done in [24], 

The evaluation of the new primed masters utilizes the methods developed in [17, 24], 
and will be explained in some detail in appendix A.2.1. The independent checks we have 
performed for the other masters, already appearing in [24], will on the other hand not be 
discussed further here, as they all produced positive results. 


3 HTL resummation 


Our goal being to extend the range of applicability of the perturbative spectral function 
to frequencies of order oj ~ gT, we will next implement in our calculation an HTL resum¬ 
mation that allows a consistent treatment of this parameter region. Following the logic of 
[17], we note that one can approximate the resummed spectral function as 

P resummed P resummed P resummed P resummed P naive P naive + Pr“Immed, (3-1) 


where ‘QCD’ refers to the full theory, while ‘HTL’ signifies a calculation performed using 
the HTL effective action. In the second step above we have used the fact that the difference 
between the QCD and HTL spectral functions should be free from IR problems, allowing 
us to perform the two calculations in a ‘naive’ form, i.e. using expansions in powers of 
the coupling g. Of the three terms on the right hand side of eq. (3.1), we have already 
computed the first one, so in the following we will only consider the naive and resummed 
HTL calculations to the necessary leading order in perturbation theory. For more details 
of the general procedure, we refer the interested reader to ref. [17]. 

The HTL calculation employs resummed gluon propagator and vertex functions, of 
which the former reads 


{Al{X)Al{Y))^^^ = 5-}Le 


abxf giQ-(X-y) 
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PluiQ) 


+ 




+ 


C Q fiQv 


g2 + nHTL(Q) Q2 + nHTL(Q) Q 


(3.2) 

where ^ denotes the gauge parameter and the transverse and longitudinal HTL self 

energies. 


Hi, 
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1 - 


Q \ ttHTL 

~T~ -i- I ’ 

iqn + qzj 


2/^2 


iQ) = 
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iqn + qz 


ml = g^N,{D - 2f [ 
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Tif 9 r-r-iO 

r 3 


- e/39 - 1 • 


(3.3) 

(3.4) 
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As discussed in section 5.2 and appendix C of [17], in the bulk channel the vertex function 
is not needed for a LO HTL calculation. In the shear case, the situation appears to be 
similar, as we have explicitly checked that the HTL vertex functions contribute to at 
least at one order higher in cu/T than the propagator corrections. We will thus ignore the 
vertex part of the calculation in what follows,^ which leads to a rather simple result for 
the resummed Euclidean correlator. 


gr(-p) 

Mac^ 


D{D-2){D-3) ^ 
2{D^-l) % 


{D + l){D-3) 


At{Q)At{Q - P) 


2(D — 2) _ - _I- 4_ _ 

^ >At{Q)At{Q-P)^ At{Q)At{Q-P) 

(^2 Q\Q-P? _o Q\Q-P? \ 

^ ^ At{Q)At{Q-P) AE{Q)AE{Q-P)i 

D{D-2){D-3) ^ f qHQ-P? _ QHQ-P? ] 

{D + 1) %\AT{Q)AEiQ - P) At[Q)Ae{Q - P) J ’ 


(3.5) 


with A^(Q) = Q^ + nHTL(Q) and A^{Q) = 


3.1 Naive calculation 

Looking first into the naive HTL calculation, we proceed to expand the expression in 
eq. (3.5) in powers of the coupling, leading us to 


gr(^) 

Mac^ 


D{D - 2){D - 3) 

L>2 _ 1 

hHTL 


Q 


{D-3){D + 1) 


Q^iQ - P? 




ng^^(Q)g^ 

Q2(Q_P)2 


(D-l) 

^ ’QHQ-Pf 

DiD-2)iD-3) ^ 

2(P-1) % 


(P-2) 


(P-3) 


QHQ-py 

q^Pn'^ 

QHQ-py 


(3.6) 


The evaluation of the sum-integrals appearing here is a rather straightforward exercise that 
utilizes the machinery developed in [17] and is explained in some detail in appendix B. The 
result of this procedure reads 


pHTL(^) 


Ma 


=- {l + 2ni-){ -h — 

dvr'' 2 M 10 45 


P 9 m 9 / 3 , 15 UJ 

- (12 H— In 2-In — 

5 V 2 2 T 


• (3.7) 


^This implies throwing out terms of order rn^uP from the result, which were kept in the bulk channel work 
of [17] and should in principle be accounted for in a consistent leading order HTL calculation. As including 
vertex corrections in the resummed HTL spectral function is, however, a technically rather tedious exercise 
and our present results indicate that even the more dominant HTL propagator corrections are negligible at 
all interesting values of cj, we have decided to refrain from performing this calculation. This issue will be 
briefly returned to in section 4. 
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3.2 Resummed calculation 


To prepare for the resummed version of the HTL calculation, we first write eq. (3.5) in the 
alternative form (cf. the discussion in section 5.4 of [17]) 


G™(P) _ D{D-2){D-3) 
Mac^ 


q^Pn^ 


,, 2 ( 112 - 1 ) ^^AT(g)AT(Q-P) 

-4 ( 72 nHTL(Q) 


- 2{D - 2) - - -4- 

^ ^Ar(g)AT(Q-P) At(Q)At(Q-P) 


(3.8) 


- (T)2 _2D-1) 


nHTL(Q)nHTL(g _ p) ^ ^^UTL^q^^HTL_ p) 


At(Q)At(Q - P) 


Ae{Q)Ae{Q - P) 


_ D{D - 2){D - 3) a; f g"ngTL(Q _ p) _ p) 

(P + l) %\At{Q)Ae{Q-P)^ At{Q)Ae{Q-P) 
where we have dropped a number of uninteresting contact terms, polynomial in the external 
momentum. Making then use of relations of the type 

1 




Qn 


A^{Q)Aj,{Q P) j 


roo j 0 

= (i + 2 n|) / — PT{q^,q)pT{^-q'^,q) 


dg° ^ ^ 0 N ^ 0 ^ '^q°'^oj-qO 


TT 


nz 


(3.9) 


where we have set P = 4 and denoted PT{qo,q) = I™{l/^T(9rn <?)}g„^.-i[go+io+] (see ap¬ 
pendix C for an exhaustive list of definitions and relations of this kind), the spectral 
function obtains the compact form {{P — Q)^ = (co — qo)^ — q^) 

32(1+ 2n^) ^ ^ f 

Jo J —o^ ^ 


4dAc2 


15(47r)2 

+ 4(?^ + Q2(4g2 ^ 7(p _ g) 

+2Q2(p _ Qfpj^{q^, q)PEi^ “ 9°; q) 

+6(g^(p - Qf + Q^{v - Q?')pTiq°,q)pEi^ - q°,q) 


q)pTiuJ - q^,q) 


V^a;-gO 

n2 


(3.10) 


All terms appearing here are of the general types already encountered in [17], and we thus 
refer the interested reader to appendix C.3 of this reference for details of their evaluation. 
The only difference between that calculation and ours is that due to the more UV divergent 
form of our integrals, we are not able to replace the distribution functions by their classical 


limits, npn^_p/n| 




in the numerical evaluation of the integrals. This implies 


that our result will explicitly depend on three scales, ca, T and m^. 

Following ref. [17], we finally write the result of the resummed HTL calculation as 

,4 




AdA 


P«tL(^) 


resummed 


AdA 


+ ^{l + 2n^)(t)l^^^{uj/T, m^/T) 


Ae 


(3.11) 


where we have for convenience separated out the naive HTL result from the rest. 
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Figure 2. The behavior of the function (j>^{uj/T) on a linear and a logarithmic scale, multiplied 


by {lo/TY in the latter case. 

4 Results 

Having now discussed the different ingredients involved in our calculation, we will next 


collect, display and analyze our results. For clarity, this section is divided into three parts: 


First, we collect the unresummed result, highlighting the fact that in the course of our 
work we discovered and corrected a subtle technical error in our original work on the same 
subject [24]. Then, we employ the HTL resummation in the form described above and 
study how this affects the spectral function, while in the final third part of the section we 


apply our results to the evaluation of the imaginary time correlator and also briefly discuss 


the shear sum rule. 

4.1 Unresummed spectral function 

As demonstrated already in [24], the unresummed NLO shear spectral function can be 
written in the form 



(4.1) 


where we have defined a dimensionless function (pVpioj/T) using the quantities intro¬ 


duced in section 2. This function was first determined already in [24], but to our surprise, 
our new formulation via eq. (2.17) was seen to lead to a different result, exhibiting a 
less divergent behavior at small oj. This discrepancy was settled after a lengthy struggle 
upon the discovery of a subtle computational error in those masters of [24] that contain a 
squared propagator. This issue, which boils down to a number of UV-suppressed analytic 
contributions having been missed in [24], is explained thoroughly in appendix A.2.2 below. 

Having settled the discrepancy, we display the behavior of the corrected (jx^{uj/T) 
function in fig. 2. A comparison with the corresponding fig. 2 of ref. [24] reveals that the 
difference between the two results vanishes quickly at large uj, yet becomes qualitatively 
important in the IR region. In particular, we see that while the analytic large-cj result 



dlTr^r® 

3a;6 


+ o(rVy 


(4.2) 
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Figure 3. The behavior of the absolute value of the (negative) shear spectral function for T = 3Tc 
(corresponding to S.TSAj^). The two blue curves stand for the NLO result evaluated with A = 
O.SAopt and A = 2Aopt, while the dashed black curve shows the leading order (LO) result. 


still holds — in accordance with the arguments of [29] — the correct w —?■ 0 limit of 
the quantity is now of order (the constant of proportionality here is not known 

analytically, nor does it have any specific physical signihcance). This has the important 
consequence that the spectral function itself has now a linear behavior at low w, implying 
that it can be straightforwardly used in the evaluation of the integrals appearing in the 
imaginary time correlator and the shear sum rule (cf. the discussion of this issue in [24]). 

Next, we insert the numerical function into eq. (4.1) in order to study the 

behavior of the spectral function numerically. Using the familiar one-loop result for cts, 
setting Nc = 3 and choosing the renormalization scale A to be varied by a factor of 2 
around the ‘EQCD value’ [37]^ 


In 


Aopt 

4 ^ 


22 ’ 


(4.3) 


we obtain the behavior shown in fig. 3. Comparing this to the results of [24], we again 
observe a fast approach of the result towards the free theory limit in the UV, but an 
important difference is that this time the spectral function does not change sign at small 
cu, but stays negative at all frequencies. 


4.2 HTL resummed spectral function 

Even though we could now directly proceed to use the corrected unresummed spectral 
function in the applications described above, it is interesting to also study the impact of 
the HTL resummation on its IR behavior, as this is expected to expand the region of 
validity of the result to frequencies of order gT. Here, our treatment is both strongly 
motivated by and closely follows the bulk computation of ref. [17], where it was seen that 
performing a one-loop HTL resummation was enough to turn the leading 0{u) IR behavior 
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Figure 4. Left: The behavior of the function (/)^rpj^(w/T, ms/T) on a logarithmic scale, with 
its value turning from positive to negative at around lo ~ l.SST. Right: The absolute value of 
the (negative) resumined shear spectral function (oj)/(A cIa) (red curves), compared with the 

unresummed NLO result already displayed in fig. 3 (blue curves) and the LO result (dashed black 
curve). For the NLO results, the two curves again correspond to the renormalization scale choices 
A = 0.5Aopt and A = 2Aopt. 


of the spectral function to an one. 

Collecting the results of the previous two sections, we see that the resummed shear 
spectral function obtains the form 


^ ^2u.y 

Ad A 2 ^ 


1 , 

10 (dvr)^ 


+ 4>%{' 




m 


iO' 




(4.4) 


where the term linear in lo in the naive HTL result (3.7) is not visible due to the way 
we chose to write the resummed HTL contribution in eq. (3.11). It, however, remains 
to be shown that the function 4)^r^-^{x,m^/T) indeed starts with an term in the IR 
limit, so that the problems we have just solved with our new unresummed calculation have 
not returned. This turns out to be a somewhat nontrivial exercise, as the function now 
depends on the extra dimensionless parameter m^/T due to the fact that we were not able 
to replace the Bose-Einstein distribution functions by their classical limits in its evaluation 
(cf. the discussion in the previous section). Inspecting the behavior of eq. (3.10) in detail, 
we see that the required linear term appears in eq. (3.11) in the limit where there is a clear 
hierarchy between the scales T and ms, ttie T. There, one can show even analytically 
that the function (j)^rj,-^{x,mE/T) indeed behaves like in the small-x limit. 

At phenomenologically interesting temperatures slightly above the critical temperature 
of the deconhnement transition, there is unfortunately no hierarchy between the scales 
T and ttie, and thus no appearance of a term linear in a; in eq. (3.10). The effects of this 
are visible in fig. d, where we display the behaviors of both /T,m^/T) and the 

entire at T = 3Tc. We observe that the HTL contribution has a visible effect on 

^As discussed in [24], at large energies it might seem more natural to choose A to be proportional to cj. 
This, however, would only have a minnscule impact on the results. 
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Figure 5. The imaginary time correlator G5j®^(r) of eq. (4.6) displayed on a logarithmic scale. Just 
as in fig. 3, altogether three functions are displayed here: Two blue NLO curves and one dashed 
red LO one. 

the result only at very small cu, but that there it again turns the result more IR divergent, 
and in fact makes the spectral function approach a constant for very small values of oj. 

As discussed already in section 3, our current resummation does not represent a com¬ 
plete LO HTL calculation, as it does not include contributions from the HTL vertex func¬ 
tion. We have, however, explicitly verified that the vertex contribution to eq. (3.7) does not 
include a term proportional to m'^ujT, but starts at earliest at order Considering 

also the fact that our present HTL result shows perfect numerical agreement with the un- 
resummed one down to frequencies a; ~ O.IT, we find it not worth the effort to undertake 
the challenge of evaluating the resummed HTL vertex contribution to the shear spectral 
function.^ In fact, due to the more bening small-w behavior of the unresummed result, we 
will choose to use it in the evaluation of the imaginary time correlator below. 

4.3 Sum rule and imaginary time correlator 

Motived by the above considerations, we will now apply the unresummed shear spectral 
function of section 4.1 to the determination of the corresponding imaginary time correlator 
and to a discussion of the shear sum rule (see e.g. [23, 27, 28]), 

1 poo q, , ^ '1 2 

- y ^ - P,{^)\t=,} = Y <T) + 0(/). (4.5) 

The linear IR behavior of our new result makes these two tasks in principle feasible, but 
the latter is still plagued by the existence of undetermined contributions to the spectral 
function of the form u!5{uj), which contribute to integrals of prj{oj)/oj. For the sake of 
curiosity, we have nevertheless performed the integral on the left-hand side of eq. (4.5) 
and found that even at order the result is within a few per cent of the energy density 

®Note, however, that for frequencies lu ~ gT, the neglected terms are still formally of 0(1) in comparison 
with terms included in (j}^{ijj/T). 
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residing on the right-hand side. This indicates that the so far unknown contact terms in 
the spectral function should be numerically subleading. 

With the imaginary time corerlator, the story is significantly simpler, as the delta 
function terms only contribute uninteresting constant terms to the integral 


G,(r) = 


r 

Jo ^ 


cosh 


— T ] UJ 


sinh ^ 


0 < r < /3 . 


(4.6) 


Plugging our unresummed spectral function to this expression, we obtain the behavior 
displayed in fig. 5. We observe an almost perfect agreement of the LO and NLO results, 
and the renormalization scale dependence of the latter is furthermore nearly invisible. 
This behavior is in practice identical to the results obtained in [24] using an IR cut-off, 
leaving all earlier comparisons with corresponding lattice and AdS results (see e.g. [36, 38]) 
unchanged. 


5 Conclusions 

The calculations reported in the paper at hand are a direct continuation of the work 
performed earlier in [24]. In this reference, we presented an NLO computation of the 
shear spectral function in pure SU(N) Yang-Mills theory, which was carried out without 
any IR resummations. While consistent with all known limits (constraining primarily its 
UV behavior), this result had one surprising feature: It suggested that the perturbative 
spectral function approaches a constant in the small-cj limit, invalidating its direct use in 
the evaluation of the imaginary time shear correlator or in sum rules. This prompted us 
to perform a further investigation of the quantity — in particular its IR behavior — by 
performing a leading order HTL resummation, expected to extend the validity of the result 
to frequencies of order gT. 

What we discovered in the course of our work was somewhat surprising: Both the naive 
and resummed HTL contributions to the spectral function turned out to produce terms 
constant in the cu —)• 0 limit, but in the presence of a hierarchy between the temperature 
and Debye mass scales, T 3> these exactly cancel each other on the right hand side of 
eq. (3.1). Puzzled by this observation, we carried out a careful and fully independent check 
of the computation first performed in [24], eventually discovering a subtle mistake there, 
described in detail in appendix A.2.2. Correcting for this was finally seen to exactly cancel 
the constant IR term from the unresummed spectral function, while leaving its earlier, 
correct UV behavior intact.^ 

Having obtained the correct unresummed and HTL resummed spectral functions, we 
proceeded to study their behavior as well as their effect on the imaginary time shear 
correlator. We observed that the HTL resummation only affects the spectral function at 
very small uj, and that the deviation of our new imaginary time correlator from the one 
derived in [24] is in practice negligible. These results can be interpreted as reflecting the 
remarkably good convergence properties of the shear channel Green’s functions, and make 

^The leading UV limit of the correlator had been predicted several years earlier in [29]. 
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US confident that we now have the behavior of the perturbative shear spectral function 
under good numerical control for a wide range of frequencies. 

With the new results at hand, the main ingredient needed before we can attempt a 
first principles extraction of the shear viscosity of SU(A^) Yang-Mills theory is obtaining 
accurate continuum extrapolated lattice data for the imaginary time correlator. We hope 
that this challenge will be tackled by several lattice groups in the near future, followed by 
coordinated lattice and perturbation theory efforts in the final analytic continuation. 
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A Master integrals in the unresummed calculation 
A.l Definitions 

The master integrals appearing in our unresummed expression for the shear spectral func¬ 
tion, eq. (2.8), are defined by 
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Following the notation of [24], we have denoted here Pt{Q) = QiiQuR^ui^) = “ (q ' P )^5 

where Pj^^{P) is the usual three-dimensionally transverse projection operator defined with 
momentum P. The sum-integration measure used here is defined as 
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•^Q go 
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d^-lq 
(27r)^-i ’ 


(A.38) 


where A and A stand for the renormalization scales in the MS and MS schemes, respectively. 


A.2 Evaluation of the masters 

The methods required for the evaluation of the above sum-integrals were largely developed 
in the bulk calculation of [17] and later generalized to the shear channel integrals in [24]. 
As most of the above cases were already considered in these two references, we will only 
discuss the new cases in detail here. In subsection A.2.1, we first evaluate the sum-integrals 
, X® , X^ and X^* that appear in the current computation due to our new way of 
treating the type (v) diagrams of fig. 1. After this, we will in subsection A.2.2 revisit 
the evaluation of a few master integrals already encountered in [24], explaining in detail a 
subtle problem we discovered in their original evaluation. 


A.2.1 New masters 

Following section 3 of [24], we write our sum-integrals in the form 



where the functions /jn read for the new masters 
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and /jo is as defined in eq. (A.37) of [24], The evaluation of these integrals follows that of 
X^, X^, and X^ step by step, except that no propagators need to be squared. We will 
thus simply list the corresponding results for the integrals below. 

-^h 

The master is related to the non-differentiated version of X^, evaluated in eqs. (B.55)“ 
(B.57) of [24]. It is easy to verify that its so-called (fz,p) and (fz,e) parts (see [24] for 
definitions) read 
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where A^- = q + {—lyEr + {—lyEqr, with Er = and Eqr = |q — r|. Finally, 

the (ps) part is obtained from eq. (B.30) of this reference by simply inserting there the 
function 


8x^ 




(A.47) 
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Following the above strategy and using eqs. (B.60)-(B.62) of [24], we get 
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In this case, the (fz,p) part can be directly read off from eq. (B.65) of [24] as 
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For the other parts, we on the other hand use the relation 
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+ ^ P, 6 ( 1 ) (^, m) + p 6 ( 2 ) (w, m) + p 6 ( 3 ) (w, m) + p @( 4 ) (w, m) , (A.52) 


and the fact that the (fz,e) and (ps) parts of the p 6 (n) (^! listed in eqs. (B.71)-(B.79) 

of [24], 
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The new master p^e^,{uj) is related to p 6 (i)(a;,m) of [24], via 
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Using this result, we immediately obtain from [24] 
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For the last of our new masters, the (fz,p) part can be read from eq. (B.82) of [24] 
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while (fz,e) and (ps) parts are available using the identity 
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A.2.2 Discrepancy with previous results 


Next, we move on to discuss the tension between our new results and those of ref. [24]. 
A lengthy independent calculation of the contributions of the different graphs of fig. 1 to 
the shear spectral function revealed a numerical discrepancy in the results of the type (v) 
diagrams. Upon closer inspection, the source of the problem was isolated to those masters 
of [24] containing a squared propagator, which indeed only appear in type (v) graphs. In 
our earlier work, our method of dealing with them consisted of introducing an auxiliary 
mass parameter in the propagator in question, and then making use of the simple 
identity 


J-=-l' d 1 

m^o I dm^ 


(A.57) 


As discussed in appendix B of [24], in the m —)• 0 limit one encounters IR divergences in 
the one- and two-dimensional integrals originating from the (fz,e) and (ps) parts of these 
masters. Our strategy with them was to separate out the divergent terms and treat them 
analytically, while the finite remainder was computed numerically after setting m to 0. 

Unfortunately, it turns out that the procedure applied in [24] missed a set of finite 
contributions that would have been correctly accounted for had we managed to set m to 
zero only after performing all the integrations. To see this in detail, consider the one¬ 
dimensional part of Pj- 3 {uj), which in [24] was written in the form 
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Of the terms here, we evaluated the integral in the last row analytically, while the terms on 
the first two rows were treated numerically after first setting m —)• 0 inside the integrand. 
The problem with this are the IR divergent terms in the subtracted part that are thrown 
out due to being explicitly proportional to m?. For some of them, the actions of performing 
the integrations and taking the m —?■ 0 limit do not commute, as can be verified e.g. from 
the analytic result 
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Similar cases also appear in the two-dimensional integrals encountered in the calculation. 


- 19 ~ 














A careful treatment of all the master integrals with squared propagators reveals that 
the results of [24] must be supplemented with the extra contributions 
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where Li stands for the polylogarithmic function. Once all of these terms are taken into 
account, one can straightforwardly verify that the discrepancy between the results of [24] 
and our new computation has completely vanished. 


B Master integrals in the naive HTL calcnlation 

The master integrals appearing in the naive HTL calculation of section 3.1 can be eval¬ 
uated using techniques familiar from [17]. First carrying out the Matsubara sums with 
standard methods and then taking the imaginary parts of the results, we obtain the simple 
intermediate results 
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where m is a non-negative integer, A a regulatory mass parameter, and Eq = ■\/q‘^ + X^. 
Taking then a derivative with respect to A (and making sure that no issues such as those 
explained in the previous section arise), we further obtain the identities 
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which leads us to the following results for the master integrals appearing in eq. (3.6): 
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As the HTL resummation is intended to be carried out in the region w <C T, we can 
simplify the above results by carrying out an expansion in positive powers of uj. For the 
integrals with the range dg, this is easily accomplished by taking the classical limit 
{\+)nq —)• T/g, which leads to the simple results 
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For the g integrals extending to infinity, we on the other hand split the integration range 
to two parts, /°° dg = / dg + dg, where A is assumed to satisfy a; <C A <C T. On 
the first of these ranges, we can again apply the ‘classical’ replacement of (l+)nq —)• T/g, 
while on the latter one, we expand the integrand in powers of uj. Taking the last term of 
eq. (B.7) as an example, we obtain for the two integration regions 
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the sum of which is clearly independent of A. This quickly leads to the results 
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which ultimately provides us with the final result for the naive HTL contribution to the 
shear spectral function, eq. (3.7). 

C Resummed HTL calculation 

In the evaluation of the resummed HTL contribution to the shear spectral function, one 
encounters the functions (cf. appendix C of [17]) 
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where we have introduced the notation 


l??l < 1, 

( 

l??l > 1, 


(C.2) 

(C.3) 

(C.4) 


Q=(g^q), Q^^(qy-q^, r/= 


(C.5) 


as well as denoted 


St(Q) = -Q^ + Ar 


mz 


\2 I ^(1 - 
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1 + 7? ’ 

[ 2 

1 — r/ 


r^(r7) = AAAI 


(C.6) 

(C.7) 

(C.8) 

(C.9) 


Upon inserting these expressions into the integral of eq. (3.10), the calculation reduces 
to one almost identical to that performed in Appendix C of [17]. We thus refrain from 
providing further details of this rather straightforward exercise, and simply display its 
numerical outcome in section 4.2. 
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